Lower limit on the neutralino mass in the general MSSM. 
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Abstract 

We discuss constraints on SUSY models with non-unified gaugino masses and Rp conservation. 
We derive a lower bound on the neutralino mass combining the direct limits from LEP, the 
indirect limits from [g— 2)«, b — > sj, B s — > u + u~ and the relic density constraint from WMAP. The 
lightest neutralino m^o ~ 6GeV is found in models with a light pseudoscalar with Ma < 200GeV 
and a large value for tan/3. Models with heavy pseudoscalars lead to m^o > 18(29)GeV for 
tan/3 = 50(10). We show that even a very conservative bound from the muon anomalous 
magnetic moment can increase the lower bound on the neutralino mass in models with [i < 
and/or large values of tan/3. We then examine the potential of the Tevatron and the direct 
detection experiments to probe the SUSY models with the lightest neutralinos allowed in the 
context of light pseudoscalars with high tan (3. We also examine the potential of an e + e~ collider 
of 500GeV to produce SUSY particles in all models with neutralinos lighter than the W . In 
contrast to the mSUGRA models, observation of at least one sparticle is not always guaranteed. 
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1 Introduction 



The minimal supersymmetric standard model (MSSM) is certainly one of the most attrac- 
tive extensions of the standard model. Apart from providing a solution to the hierarchy 
problem it also fits in very well with the idea of the gauge unification and gives as a bonus 
a good dark matter candidate represented by its lightest supersymmetric particle, the 
LSP for short. One the other hand, to the question of how light this LSP can be so that 
one has an idea about the masses for the rest of the spectrum, the model provides no an- 
swer. The reason for this is that there is no definite model for the masses of the sparticles 
(and other scales) which would ultimately point out to a theory of the (soft) breaking of 
super symmetry. One exception relates to the mass of the Higgs if one makes the mild 
requirement that the scale of supersymmetry breaking to be much above the electroweak 
scale. This gives one of the robust predictions of the model which is that there must be 
at least one light Higgs with mass below about 130GeV [JJ, for alternatives see [21 E]. For 
practically all the other (s)particles of the model one does now have some lower bounds 
on their masses from extensive phenomenological studies and searches both at colliders 
and in astroparticle experiments. One can for example mention the impact of both LEP 
and the Tevatron in constraining the model parameters through the direct searches, the 
importance of rare decays notably b — > S7 or even low energy experiments like the recent 
measurement of the muon anomalous moment (g — 2) M Cosmological considerations 
like the thermal relic density of the LSP also put some stringent bounds on the MSSM. It 
should however be emphasised that many of the bounds and constraints, especially those 
from indirect searches, but not only, are based on theoretical assumptions some of which 
with strong prejudices. The aim of this paper is to critically review the lowest bound on 
the neutralino LSP with as little theoretical bias as possible. Instead of reviewing a host 
of models, we will try to delimit the parameter space that has an incidence on the LSP 
mass bound in a, as much as possible, model independent way 0. We will see that one 
does not have to play with a large number of parameters to arrive at the lowest bound 
on the LSP, instead the parameter space will be naturally reduced to only a few critical 
parameters. 

A large amount of work has been devoted to the so-called minimal supergravity model 
(mSUGRA) mainly because of its simplicity. Indeed, instead of dealing with about 100 
parameters in the general MSSM, one works with only four parameters and a sign. With 
such a constrained model it becomes possible to combine the various experimental data 
and severely delimit the parameter space. Within this model the lowest limit on the LSP 
is extracted in a straightforward way from the LEP data which give m^o > 59GeV [01 Ej . 
Moreover, with the new precise determination of the relic density of cold dark matter by 
WMAP [Hj, the mSUGRA parameter space considerably shrinks to thin strips [Hj suggesting 
a rather focused phenomenology and benchmarks for the next colliders [TU| . Supersym- 
metric signals can be radically different in other scenarios, scenarios where one can have 
much more room for manoeuvre even after imposing all the latest data. In particular, the 
LSP bound can be much smaller than the one found in the context of mSUGRA as we 
will see. 

Different approaches have been followed aiming at a reduction of parameters. Theoreti- 
cally motivated approaches include the so-called anomaly mediated (AMSB) models ^JJ 
and generalisations J2j for example, or schemes that appeal to superstrings that are ei- 
ther moduli-dominated or contain a mixture of moduli and dilaton fields fSl EH E3 ITT)] 
with varying degrees of simplification. There are even scenarios that interpolate between 
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mSUGRA and some of the AMSB and string inspired methodology ^JJ. Another ap- 
proach without necessarily an underlying framework for supersymmetry breaking, lifts, 
at the GUT scale, the mass degeneracies of the mSUGRA model either in the sector of 
the gaugino [HI CHI E01 EU E2] , Higgs [23 El] or the sfermions (2HI ETJ . Some analyses let 
loose the shackles of mSUGRA and models defined at high scale but still impose ad-hoc 
constraints, like equality of scalar masses at the electroweak scale [26 . Reviewing this 
situation one can still ask about the lowest stable LSP mass allowed by the present data 
regardless of a pre-defined scheme. This is important because not only one would want to 
know what the impact of the lightest LSP is for cosmology and direct detection but how 
the collider phenomenology can get affected when one puts aside assumptions related to 
the LSP mass. For example one might inquire about the sensitivities of the direct search 
detectors such as CDMS [27], EDELWEISS j2H] and ZEPLIN [25] for neutralino masses below 
those typical of the LSP of mSUGRA. The LSP, in the R-parity conserving MSSM, being 
the end product of the decay of any super symmetric particle, is important for signatures 
at the colliders since it will always show up in any analysis. Another reason why it is 
interesting to study the light neutralino is that it opens the door for a sizable branching 
fraction of the Higgs into invisible thus reducing all the branching fractions into the usual 
discovery channels, h — > 77 or h —>■ bb |301 I3T] . This possibility has triggered analyses 
by the LHC experiments to search for the Higgs with a sizable branching fraction into 
invisible [32*] 133] . This has also motivated studies to reanalyse Tevatron data in models 
with light neutralinos and establish whether or not supersymmetry could be discovered 
there in the chargino/neutralino channel leading to tri-leptons 

The oft-quoted lower limit, m^o > 59GeV [HIE] which applies to mSUGRA, is basically 
derived from the quite robust model independent lower limit on the chargino mass m-+. 
The (lightest) chargino mass is set from chargino pair production at LEP2. Deriving the 
lower limit on the neutralino mass from this analysis, on the other hand, tacitly assumes 
a model of unified gaugino masses at the GUT scale, much like what is hypothesised in 
mSUGRA. However in a general MSSM, the masses m f o and m~+ are uncorrelated and 
the lower limit on the LSP neutralino mass weakens considerably if only the LEP constraint 
is taken into account. This occurs for a bino neutralino when Mi <C Mi where M\(M2) is 
the U{1)(SU{2)) soft-susy breaking gaugino mass. To derive a lower bound on the LSP 
in these conditions one has to turn to cosmology. Indeed, the neutralino LSP cannot be 
too light as it would conflict with the precise measurement of the thermal relic density of 
cold dark matter inferred from WMAP [8J. One is then indirectly led to include the sleptons 
into the picture, in particular the right handed ones, Ir. Indeed a light neutralino that 
is mostly a bino will annihilate preferably into fermions through right-handed sleptons 
because the latter have the largest hypercharge among all sfermions [3J3 EH] . As a rule 
of thumb, with all sfermions heavy but the three right sleptons, a rough approximate 
requirement is 



with all masses expressed in GeV. We identify £7 X with the fraction of the critical energy 
density provided by neutralino LSP and h is the Hubble constant in units of 100 km sec -1 
Mpc -1 . Therefore the lower bound on the slepton from LEP2, when used in conjunction 
with the upper limit from the relic density 1 , plays an important role in constraining a 
light neutralino. We will see that to accommodate the relic density constraint, a slight 

1 Note that we only consider scenarios where the neutralino contributes to the cold dark matter. 




(1) 



2 



Higgsino component can make the mostly bino LSP couple to the Z and the Higgses. If 
one is not far from these resonances there will be efficient annihilation of the neutralino 
LSP. 

To further constrain the lower bound on the LSP mass and bring in new observables in 
the picture, some assumptions on the sfermion masses need to be made. A well motivated 
mild assumption relies on so-called minimal flavour violation in order to evade bounds on 
flavour changing neutral currents. In the context of a generic minimal flavour violation 
model, sfermions with the same quantum numbers would share the same mass, with per- 
haps a slight breaking from Yukawa mixing. We would then treat all the sleptons on an 
equal footing. This then brings into the picture the constraint from the measurement of 
the muon anomalous magnetic moment, 5a fl . The latter is also sensitive to the presence 
of light charginos/neutralinos and smuons. Both the theoretical predictions and exper- 
imental results on the muon anomalous magnetic moment have been refined on several 
occasions in the last year and the situation is still evolving. We will therefore be very 
conservative in imposing limits from (g — 2)^ using numbers that encompass different 
estimates. Nevertheless we will find that the upper bound on the (g — 2)^ does constrain 
some of our scenarios. In our study we will also show how results change if one does not 
impose the (g — 2) M bound. 

The squarks on the other hand hardly enter the game. One needs however to assume 
them to be sufficiently heavy so that, especially for large tan/3, one does not conflict 
with the bound from b s'j. Moreover, squarks of the third generation should be heavy 
enough and or the tri-linear mixing parameter of the top be large so that, especially 
for low tan/5, one evades more easily the lower bound on the Higgs mass. One usually 
also chooses the pseudoscalar Higgs mass to be heavy enough to achieve this. With the 
assumption that the pseudoscalar Higgs is heavy (beyond say ~ 300GeV), an absolute 
lower limit on the neutralino mass, m^o = 18GeV(29GeV) for tan/5 = 50 (10) results 
from combining cosmological and collider constraints. 

Relaxing the assumption that the pseudoscalar is very heavy, after all the LEP lower 
bound is only about 90GeV, leads to an even smaller lower bound on the neutralino mass, 
m^o > 6GeV. This is because, as was pointed out in [3B], a new possibility for escaping the 
relic density constraint for very light neutralinos opens up: neutralino annihilation into bb 
pairs via a Higgs resonance. This channel becomes efficient in the large tan (3 regime due to 
the enhanced couplings to b quarks and r's. We will show that the concomitant presence 
of a light charged Higgs at large tan/3 means that these models are tightly constrained 
by b — > 57 as well as by B s — > n + fi~. Nevertheless there remains a possibility for very 
light neutralinos, m^o ~ 6GeV, when tan/3 > 25. The phenomenological implications of 
these classes of models is markedly different from the models with heavy pseudoscalars 
and deserves special attention. Several Higgses with masses in the 100 — 200GeV range 
and with high tan/? can be detected at the Tevatron , even with 2fb~ 1 . These models 
are testable at the current Runlla. In this scenario, the Higgs exchange also dominates 
the spin-independent neutralino-nucleon cross-section. Large cross-sections are found for 
very light neutralinos, m^o ~ 6 — lOGeV. 

This paper is organised as follows. We first critically review the constraint from LEP, 
paying attention both to the limits on neutralinos, charginos as well as sleptons. We 

Unstable neutralinos with masses in the MeV range have been entertained |37j and their astrophysical 
(supernova) implications studied. 
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will then analyse the constraints from the relic density and from the muon anomalous 
magnetic moment, (g — 2)^ and the rare decay b — > S7. Other constraints that we use, 
although in most cases without much impact, are the effect on the invisible width of the Z, 
Z — > bb and the rare decay B s — > The latter can play a role for a light pseudoscalar. 

We will concentrate on the intermediate to large tan f3 region and present absolute lower 
limits on the mass of the neutralino LSP. The case of models with a light pseudoscalar 
Higgs deserves a special section both as concerns the impact of the various constraints 
as well as the prospects at future colliders. Here the Tevatron plays an important role. 
The prospects for direct detection experiments will be discussed in this context. Finally, 
a section will be devoted to the potential of the linear collider to produce supersymmetric 
particles in models with light neutralinos, m-+ < M w . We will close with a conclusion 
that summarises our results. 



2 MSSM parameters 
2.1 Physical parameters 

When discussing the physics of charginos and neutralinos it is best to start by defining 
one's notations and conventions. All our parameters, unless stated otherwise, are defined 
at the electroweak scale. The chargino mass matrix in the gaugino-Higgsino basis is 
defined as 



M 2 V2M W cos (3 

y/2M w sin (3 ji 



(2) 



where M 2 is the soft SUSY breaking mass term for the SU(2) gaugino while fi is the 
so-called Higgsino mass parameter whereas tan f3 is the ratio of the vacuum expectation 
values for the up and down Higgs fields. 



Likewise the neutralino mass matrix is defined as 

/Mi -MzsinfVcos^ M z sin 9 W sin (3 \ 

M 2 M z cos 9 W cos P —M z cos 6w sin (3 

— M z sin9 w cos/5 M z cos 6w cos (3 — /x 

\ M z sin 9w sin P — M z cos^iy sin/3 —\i / 



(3) 



where the first entry Mi (corresponding to the bino component) is the U(l) gaugino mass. 
The oft-used gaugino mass unification condition corresponds to the assumption 

Mi = - tan 2 9 W M 2 ~ M 2 /2 (4) 
3 

Then constraints from the charginos alone can be easily translated into constraints on 
the neutralino sector. Relaxing Eq. HI or removing any relation between Mi and M 2 
means that one needs further observables specific to the neutralino sector. It is a trivial 
observation that if Mi <C M 2 , one can get a very low neutralino mass independently 
of the chargino mass derived from Eq. |21 

In the approach we are taking, the free parameters include the ones of the gaugino 
sector as well as the parameters of the slepton sector. As discussed earlier, we will 
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assume that all the squarks are heavy. Allowing light squarks would not affect much the 
cosmological constraints for a neutralino that is mostly bino and in any case light squarks 
are strongly constrained from b — > sj. For the slepton sector we will consider primarily 
a model reminiscent of mSUGRA models but without gaugino mass unification. This 
model features a common mass for the sleptons at the GUT scale. As we will see, for the 
purpose of constraining the LSP mass, the slepton mass spectrum is essentially the same 
had we imposed the more general requirement of minimal flavour violation by giving a 
common mass to the left and a common mass to right sleptons of all three generations. 

Using the renormalisation group equations, the masses at the weak scale can be related 
to the ones at the GUT scale. For sleptons this can be done rather independently of the 
other MSSM parameters. Defining the parameter 

which characterises the amount of non universality in the gaugino masses, the weak scale 
slepton masses write 

ml R = Ml + .88 r 2 2 M 2 - sin 2 9 W D Z 

ml L = Mq + (0.72 + .22 r 2 2 )M 2 - (.5 - sin 2 B W )D Z 

ml e = Mq + (0.72 + .22 r 2 2 )M 2 + D z /2 with 

D z = Mfcos(2/3) (6) 

Here the gaugino masses, M 12 are defined at the electroweak scale. To obtain M 1 = r u M 2 
at the electroweak scale one needs Mi ~ 2r 12 M 2 where M, are defined at the GUT scale 
and M 2 « 0.825M 2 . In the mSUGRA models one has Mi = M 2 . The parameter M is 
here defined at the GUT scale. For the sake of simplicity we have neglected the Yukawa 
couplings in the RGE that could affect the f sector, however mixing will be taken into 
account through the \x term as we will see. 

We will be mainly interested in models where 7*12 << 1, then in the RGE equation all terms 
in 7*12 will be negligible. One obtains a natural splitting between the right-handed/left- 
handed sfermion masses. Indeed, m e - fl w Mo and is typically much lighter than cl which 
receives in addition a contribution, m§ oc M|. 

The non universality in the GUT-scale relation for gauginos which we investigate 
in this paper are quite plausible as many models beyond mSUGRA feature some kind 
of non-universal masses. For example SUGRA models with non-minimal kinetic terms 
[TH] , superstring models with moduli-dominated or a mixture of moduli and dilaton fields 
[T5J H3J H2J [13 HZ] and anomaly mediated SUSY breaking models [TT] all feature non- 
universal masses in the gaugino and/or scalar masses. 

Note in passing that Eq. |U] can be extended to squarks and if we take M 3 = r 32 M 2 with 
r 32 > 1 at the GUT scale one could make the squarks "naturally heavy" as we have 
assumed. Then the gluino mass relation m g ~ 4M 2 obtained in mSUGRA type models 
turns into rrig ~ 4r 32 M 2 (r 32 > 1). First and second generation squark masses, neglecting 
the small bino contribution oc M 2 , can be approximated as 

"4,* ~ m Ln + °- 6m l ( 7 ) 

We fill first consider the case where the pseudoscalar is heavy, Ma = 1 TeV. Altogether 
we allow 6 free parameters at the electroweak scale and fix M 3 = rrig = 1 TeV : 

tan/?, M x , M 2 , M , /1, A, (8) 
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A q j are the tri-linear terms for the quark q and lepton I. For scans over the parameter 
space, unless otherwise specified, we will consider the range 

5 < tan/3 < 50, 
M 2 < 2000 GeV, 

.001 < r i2 < .6, 
< 1000 GeV, 
\A t \< 2400 GeV, 

M < 1000 GeV (9) 

We will fix Ai = for all sleptons as most of the processes we will discuss are not very 
sensitive to the exact value of this parameter. Although the mixing in the stau sector, 
which is oc A T = A T — /itan/3, can be relevant for the calculation of the relic density, 
in our case we always have \fi\ tan/3 > 1000 so this term usually dominates the mixing. 
Finally we will consider also the case where the pseudoscalar mass is a free parameter 
varying it in the range 

92GeV <M A < ITeV (10) 



3 Limits from LEP 

The direct limits from LEP2 on charginos, neutralinos as well as on sleptons are relevant 
for the lower bound on the lightest neutralino. As argued, the sleptons play an important 
role in the relic density calculation especially when the light neutralino is a bino. Here we 
revisit the lower limits that can be obtained by LEP2 on sfermions and gauginos/Higgsinos 
when one relaxes the gaugino universality assumptions. 

Charginos 

The lower bound on the chargino mass rests near the kinematic limit, m+ > 103.5GeV 
unless the sneutrino lies in the range 75 < m fe < 85GeV, then the bound drops to 
m+ > 73 GeV 7\. This is due to the destructive interference between the t and s- 
channel contributions to e + e~ — > XiXi ■ These bounds basically translate into bounds on 
(M 2 , > (73,103)GeV. 

Neutralinos 

The LEP experiments quote a lower limit on the neutralino mass, m^o < 59.QGeV, 
while assuming unified gaugino masses at the GUT scale, r i2 = 0.5 This constraint 
on the neutralino mass is basically derived from the lower limit on the chargino mass 
obtained in the pair production process, e + e~ — > XiXi which depends on M 2 ,/i. It is 
only through the relation between the gaugino parameters Mi and M 2 that the neutralino 
mass limit is obtained. In a general MSSM, the charginos and neutralino masses are 
uncorrelated and the lower limit on the neutralino mass weakens when r 12 < .5 becoming 
roughly m^o > r 12 x m^+. Taking into account the constraint from the chargino sector, 
the lightest neutralinos are then typically found in scenarios where Mi <C M 2)/ u, that is 
the light LSP mass is set by the parameter Mi and is thus mostly a bino. 

In such scenarios, the processes e + e~ — > Xi/dL X1X3 can De used to somewhat constrain 
the parameter space. In particular, when sleptons are light, the neutralino cross-sections 
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depends crucially on the Higgsino content. The cross-sections are generally enhanced and 
regions with small \i can be more severely constrained by the neutralino production cross 
sections than by the chargino process [30J. In our scans, we have implemented the upper 
limit from the L3 experiment on a(e + e~ — > X1X2 + X1X3 ~^ with I = e,/i, using 

the tables in jSH] that simulate both the signal and background. 

The radiative processes where a photon is emitted in addition to a pair of invisi- 
ble supersymmetric particles, like the lightest neutralino, will contribute to the process 
e + e~ — > 7 + invisible which has been searched for by the LEP2 experiments. The neu- 
tralino radiative process e + e~ — > X1X1I1 is essentially dominated by the t-channel selec- 
tron exchange diagrams and has a small cross-section El| • We computed the double 
differential cross-section(energy and angular distribution of the photon) as well as the to- 
tal cross-section exactly using calcHEP [UJ for a wide range of parameters in the MSSM. 
The cross-section can reach a = 50 fb for m e - B « lOOGeV but drops steadily as the scalar 
mass increases. While a few events are expected at LEP2, this is not enough to overcome 
the uncertainty in the standard model neutrino contribution to the 7+invisible channel. 
Only mild improvement is obtained from fitting to the full energy-angular distribution as 
the distribution is somewhat similar to the background, peaked near low-energy photon. 
We conclude that one cannot constrain any further the neutralino mass using this process 
independently of the selectron mass. 

In principle the upper limit on the invisible width of the Z, Y 'z inv < 3MeV from 
LEP1 can also constrain the lightest neutralino. In the particular case of a bino LSP 
with Mi <C M2, the coupling to the Z strongly depends on the amount of Higgsino 
mixing. Thus one obtains a lower bound on /1, in the large tan/? limit, > HOGeV. 
This constraint is in general already satisfied after taking into account the LEP2 data on 
chargino/neutralino production. 

Sfermions 

For selectrons, a limit of 99.5GeV can be set on both e~i as well as cr in the case 
of a light neutralino, whereas basically model independent limits of > 96GeV and 
TOf > 86GeV can be reached [Zj. It is important to point out that the lower limit on the 
stau is about lOGeV smaller compared to the smuon and the selectron as this will have 
a consequence on the relic density contribution. Note also that because of the mixing is 
stau sector (which in our case is induced solely through the \i term), we can with the 
same minimum value of M in Eq. El arrive at lower values for staus than for selectrons. 

Higgs 

In models where the pseudoscalar mass is heavy, the limit on the lightest CP-even 
Higgs mass from LEP2, M h > 114.4 GeV, applies. We have used FeynHiggsf ast [12] 
to calculate the Higgs mass and have imposed the limit Mh > 113 GeV to allow for 
theoretical uncertainties. 

In models with a light pseudoscalar the above LEP2 constraint is relaxed. When 
Mh ~ Ma and cos(a — 3) ~ 1, the channel e + e~ — ► hZ is strongly suppressed. LEP2 can 
only make use of the hA —>■ bbbb, rrbb channels to put an absolute bound of Mh, Ma > 
91.6 GeV [4*3] 144[ 143] . In these models, the heavy CP-even Higgs channel e + e~ — ► HZ 
which is oc sin 2 (a — (3) is favoured. Unfortunately the mass splitting between the two 
scalar Higgses (Mh > Ma) is generally sufficient to put the heavy Higgs beyond the reach 
ofLEP2. 
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4 Indirect limits : Q/i 2 , (g — 2)^, b — > sj, B s —> 
and Z —> bb 

4.1 Relic density of neutralinos 

The MSSM model with a light stable neutralino must be consistent with at least the 
upper limit on the amount of cold dark matter. Here we take the latest bound from WMAP, 
tth 2 < .128 and also compare with the old limit (of 2001) from BOOMERANG 01,, MAXIMA 
|Hj and DASI [IB] with Vth 2 < 0.3. We will refer to this limit as pre-WMAP. Our calculations 
of the relic density is based on micrOMEGAs 49 , a program that calculates the relic density 
in the MSSM including all possible coannihilation channels jlHj. For the light neutralino 
masses under consideration, it is the main neutralino annihilation channels that are most 
relevant, in particular annihilation into a pair of light fermions. Basically two types of 
diagrams contribute, s-channel Z (or Higgs) and t-channel sfermion exchange. A light 
neutralino that is mainly a bino couples preferentially to right-handed sleptons, the ones 
that have the largest hypercharge. To have a large enough annihilation rate (in order 
to bring down the relic density below the upper limit allowed) one needs either a light 
slepton or a mass close to Mz,h,H,A/2- In the former case, the constraint from LEP on 
sleptons and in particular staus plays an important role. In the heavy slepton case, the 
coupling of the Z should be substantial, which requires that the neutralino should have 
a certain Higgsino component |3lH I31j. This means \x small, but still consistent with the 
chargino constraint. 



4.2 (g - 2), 

To derive the bound on 5a a = af, xp ' — a* heo we take into account the recent correction 
to the sign of the hadronic light by light contribution, a^ BL |50[ I51j. Considering that 
there are still a few issues that need to be clarified, see for instance [H2], we allow a larger 
theoretical error and take like in [S3] a^ BL = (+80±40) x 10 -11 . For the hadronic vacuum 
polarisation we take an average of the new results. Jegerlehner has reported a new value 
that includes BES, CMD-2 with a value on the hadronic part a^ ad = (6889 ± 58) x 10~ n 
[5%] . Teubner [S3] quotes two values (inclusive and exclusive, depending how the data is 
integrated) but uses QCD sum rules to favour the value extracted from inclusive data, 
a had _ gg3i -j- 5g _|- 20) x 10~ n . Finally Davier reported a value based on e + e~data 
a^ ad = (6847 ± 60 exp ± 36^) x 10~ n as well as another one based on r + r~ data, a^ ad = 
(7090 ± 47 ± I2 exp ± 38sf/(2)) x 10~ n . Since these two values are slightly inconsistent we 
will first consider a limit coming from an average of the hadronic estimates from e + e~ 
alone as well as a more conservative limit encompassing both e + e~ and r results. 2 

The theoretical estimate is obtained after adding the pure QED contribution, a® ED , 
the weak contribution, a™ eak and the three different hadronic contributions: a^ ad ,a^ BL 
and the NLO hadronic contribution, a^ ad//NLO , 

theo. _ QED , weak , had , LBL , had/NLO 



2 The new update analysis of Davier et al..|56| finds a similar r-based result and a much better agree- 
ment between the e + e _ -based and r-based value. We have not used these values in our analysis as our 
conservative range encompasses already the newest estimates. 
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With the latest experimental data on the g — 2 measurement [2. , bringing the world 
average to 

a ex P . = 11659 203 ± 8 x 1(T 10 (11) 
we get, after averaging of the three different hadronic contributions from e + e~ data alone, 
6a^ = a^ xp - - af eo - = (33.6 ± 8, cxp . ± 11.5| thco .) x 1(T 10 (12) 
Adding linearly the theoretical error to a 3a experimental error leads to the allowed range, 

- 2 < 5a^ x 10 10 < 69 (13) 

Including the r data analysis of Davier (HE] in the average, reduces the discrepancy 
between theory and experiment 

8a p = (27.8 ± 8, cxp . ± 11.3| thco .) x 1(T 10 (14) 

The 3cr allowed range now reads, 

- 7.5 < 5a^ x 10 10 < 63 (15) 

As the results have changed frequently in the last months and since many issues need to be 
clarified in particular in the estimation of the hadronic contribution, the most conservative 
allowed range is found by including two allowed bands, instead of performing an average. 
Combining the value obtained from e + e~ data (Eq. fH?J) with the one from r data alone 
of Davier, one gets a 3a range 

- 25 < Sa^ x 10 10 < 69 (16) 

where the lower bound is set by the r data and the upper bound from Eq. We will 
refer to Eq. E3as our conservative bound on (g — 2) M . 3 

In the MSSM, a M gets a contribution both from neutralino and chargino loops, the 
latter being dominant. One expects large effects for large tan f3 and light sfermions {Jxjv^) 
[H7j . Typically, the sign of \i is strongly correlated to the one of 5a^. However there 
are special cases where cancellations can occur between the chargino and the neutralino 
diagrams thus flipping the relative sign of 5a fl and /i. This happens for charginos much 
heavier than neutralinos and the resulting value of 5a^ is small. This property will open 
a small corner of parameter space where models with /i < will be allowed even when 
imposing the bound in Eq. ED When fi > on the other hand, one expects rather mild 
constraints on the parameters of the MSSM, furthermore the constraints will be mostly 
in the chargino/slepton sector rather than directly on the neutralino sector. 

4.3 b -> S7 

Our calculation of the b — > 37 closely follows the approach of Kagan and Neubert jSH] (see 
also 59 ) with NLO, bremstrahlung and some non-perturbative effects taken into account. 

3 This last bound is similar to the one we had obtained (at 2a) |S] before the new calculations and the 
more precise results came out in summer 2002 4 . 
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This approach makes it easier to include the effects of New Physics. Our rates refer 
to the so-called "faked" total rate with a cut on the photon energy corresponding to 
5 = 0.9 (E 7 > m b /20). We have updated this analysis by including in a simple way the 
recent suggestion of Misiak and Gambino jHO] of using the MS charm mass. The Wilson 
coefficients for the standard model and the charged Higgs are evaluated at the NLO [61, 
whereas we only include the remaining SUSY contributions at LO with however the 
important inclusion of the enhanced large tan (5 effects [HBj (SUSY threshold corrections 
to the running b mass). For the latter we include both the strong a s contribution as well 
as the Yukawa contribution for the SUSY Wilson coefficients as well as the charged Higgs 
and the Goldstone contribution 4 . Our standard model value (with scale parameters set 
at nib), gives Br(b — > sj) = 3.72 x 10~ 4 while the scale and other parameter uncertainty 
(a s , CKM matrix elements) are about 10%. To bound the SUSY contribution we take 
the world weighted average of the CLEO jHlj BELLE 65J and ALEPH [HE] measurements 

Br(b -> sj) = 3.23 ± .42 x 10" 4 (17) 

We require that after allowing for the (scale) uncertainty in the theory calculation the 
result must be within 2a of the experimental result, Eq. El Since the theory uncertainty 
is roughly constant over the SUSY parameter space and in order to have a faster scan 
we have allowed for a conservative fixed uncertainty of 10% independently of the SUSY 
parameters. Thus in effect we require the theory prediction to fall within the range 

2.04 < Br{b -> s 7 ) x 10" 4 < 4.42 (18) 

In the MSSM, contributions to the b — > s 7 depends mostly on the squark and gaug- 
ino/Higgsino sector as well as on the charged Higgs. The heavy squarks that we consider 
do not completely decouple and one can get substantial corrections to the SM branch- 
ing ratio. In particular at large tan f3 there is a strong A t dependence through a term 
oc At tan (3 from the mixing in the stop sector. However, the light Higgs mass is also sensi- 
tive to the mixing in the stop sector. As long as the pseudoscalar is heavy, one finds that 
for small values of tan/3, one can find values of A t , typically A t large and positive, that 
satisfy both the Higgs mass limit as well as the b — > s 7 bound. For light pseudoscalars 
on the other hand we will see that it is very difficult to satisfy the b — > 57 bound. At 
large values of tan f3, the Higgs bound is more easily satisfied and one can pick a set of 
values for A t (typically A t < 0) that allows the correct amount of b — > s 7 whether or not 
the pseudoscalar is heavy. Note that the light squarks contribution to b — > 57 rapidly 
becomes too important. This again justifies our choice of large squark masses. 

In the end we find that because of the free mixing parameter in the stop sector, the 
b — > s 7 constraint has a mild impact on the models considered at least as concerns the 
mass of the neutralino. One exception is the light pseudoscalar case that will be discussed 
in more details in the last section. 

4.4 B s — > 

The CDF experiment at Fermilab has obtained an upper bound on the branching ratio 
B.R.(B S -> < 2.6 x 10" 6 ) [HH and should be able to reach B.R.(B S -> n+fj,- < 

4 We have corrected some typos in . We thank Paolo Gambino for checking our results and agreeing 
with our implementation. 
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2. x 10~ 7 ) in Runlla. In the SM, this branching ratio is expected to be very small 
(~ 3 x 1CT 9 ). In the MSSM SUSY loop contributions due to chargino, sneutrino, stop 
and Higgs exchange can significantly increase this branching ratio. In particular, the 
amplitude for Higgs mediated decays goes as tan/5 3 and orders of magnitude increase 
above the SM value are expected for large tan/3. This process is relevant mainly in the 
light pseudoscalar scenario. Our calculation is based on |HH] and agrees with [20] • Amj 
effect relevant for high tan/? are taken into account. 

4.5 Z -> bb 

We have also included the constraint from Z — > bb although in most cases it is completely 
harmless. For example, the charged Higgs also contributes to Z — > bb. However this 
puts a constraint on the Higgs sector only for Higgs masses below the LEP limits. Our 
calculation of Z — > bb is along the lines of [7U] but we have corrected a few typos contained 
in [70]. 

5 The lower bound on m~o when Ma = 1 TeV 

Al 

In this section we present our results for the case of a heavy pseudoscalar after imposing all 
direct and indirect constraints. In particular we have imposed only the conservative limit 
on (g — 2) M , Eq. ED Nevertheless it is worth discussing the impact of the much stricter 
bound, Eq. E3 As the allowed region in \i < models changes drastically depending on 
how one evaluates the hadronic contribution, we will defer the extensive discussion of this 
constraint in the section devoted to the \i < case. 

5.1 fi>0 

We first consider the case /i > 0. In these scenarios, the relic density of dark matter 
provides the main constraint on models with light neutralinos. Indeed the LEP limits on 
charginos and neutralinos allow very low values for neutralino masses provided M% << M 2 . 
The (g — 2) M constraint is in general easily satisfied except for very large values of tan f3. 
The constraint for b — > 37 also affects mostly models with large values of tan (5. 

1) Heavy sleptons: M = 500GeV. 

In this situation all sfermions are heavy. As a consequence their contribution to the 
relic density is negligible. Annihilation of the LSP dark matter through s-channel Z 
and Higgs exchange is on the other hand possible. Then in order to have a sufficient 
neutralino annihilation rate, to satisfy the upper limit on the relic density, a large enough 
coupling to the Z is necessary. This requires the LSP (which is mostly a bino) to have a 
certain amount of Higgsino component, especially as one moves away from the Z peak. 
Bino-Higgsino mixing has little tan/5 dependence and scales like Mz/n which calls for 
the smallest \y\ possible. It can also be enhanced somewhat if there is little splitting 
between M\ and fi, but this would generally not minimise the LSP mass. We scanned 
over the parameters ri 2 ,M 2 ,/i and tan/3 as specified in Eq. |H1 The relic density as a 
function of the neutralino mass clearly shows the effect of the Z peak and the lower limit 
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Figure 1: a) Relic density of the LSP vs m^o in the MSSM with Mq = 500 GeV. To guide 
the eye the WMAP upper bound is displayed, b ) Values of \i consistent with the LEP and 
relic density constraints. We scan on tan/3,M 2 ,r 12 and [i. 



on the neutralino mass, m^o »s 32Ge^, Fig. [T^.. Note that although some points give a 
relic density that is much too low, models with Qh 2 w 0.1 can be found for any value 
of the neutralino mass consistent with the lower limit of ~ 32GeV. In particular for the 
lightest neutralinos in this scenario (say within 2GeV of the lower limit) the value of the 
relic density falls within the range of WMAP. Incidentally, Fig. ^also shows that the relic 
density can drop quite dramatically even beyond the Z-peak and upon inspection even 
higher \x values are allowed in this range. This corresponds to the contribution to the 
s-channel lightest Higgs h with m^o pa m^/2. Since in the scenario with Ma = ITeV, 
rrih > 113 GeV this occurs beyond the Z-peak. One of the reasons this second peak looks 
broader is because we are scanning over a wide range of h masses. Another reason is that 
the coupling of the (almost) bino LSP is larger for the Higgs than it is for the Z. For the 
latter the coupling is in fact quadratic in the Higgsino-bino mixing. As we will discuss 
later, we expect the presence of light Higgses other than h to reduce the value of flh 2 and 
for masses of the order of the Z mass to allow an even lower limit of neutralino LSP. 
Fig. clearly shows the preference for a significant Higgsino component, the small [i 
region, when the neutralino mass is below Mz/2. The minimum value for the LSP mass, 
m^o p^ 32GeV, occurs for ri 2 < .2, see Fig. El As rvz approaches the mSUGRA value, 
ri2 = 0.5, the limit from the chargino mass at LEP dominates and the relic density 
constraint has no effect. This lower bound is more or less independent of tan/3 in the 
range under consideration 5 < tan/5 < 50. Fig. El also shows how the improvement of the 
bound on the relic density from WMAP has strengthened the lower limit on m^o by pa 5GeV 
when r 12 < .2. For larger values of the non universality parameter, the improvement 
from WMAP is only marginal. Indeed, already with flh 2 < 0.3, one is confined to a region 
(m^o > 35 GeV ) close enough to the Z peak so that flh 2 drops sharply with m^o (Fig.[T]a). 

The preference for the small \x region also implies that light neutralinos, say below 
40GeV, are necessarily accompanied by light charginos (below 250GeV). However, as soon 
as one moves closer to Mz/2, the Higgsino component does not have to be too small and 
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Figure 2: Lower limit on m^o vs r 12 after imposing the relic density constraints. The 
limit from LEP2 alone is also displayed. Here M = 500 GeV, \i > and we scan on 
M2,tan/3, \fi\. The region below the lines is excluded. 

one generally has m-+ < 450GeV. This has important consequences for colliders. For 
example it could be difficult for the Tevatron to find charginos in the trilepton channel 
{x + xV)- Similarly a 500GeV collider would not be able to find charginos in some scenarios 
as will be discussed in the last section. Note that under the conditions specified here, no 
improvement on the lower limit on the neutralino can be set from (g — 2) M . 

2) Light sleptons: M < 500GeV 

In addition to the effect of the Z exchange one expects the contribution from t-channel 
sfermions to XX ff to weaken the constraint from the relic density. Here the constraints 
on the masses of selectrons/smuons and staus from LEP2 are important. The latter is 
relevant for large values of tan (3 when the f can be much lighter than the selectron. As 
the left-right mixing in the stau sector is proportional to A T — /i tan (3 this is especially 
true in the large fi region. The (lightest) f contribution to the relic density in setting 
the lower bound on the LSP is largest among all sleptons. Not only there is a S"-wave 
contribution (oc m 2 ) but the (dominant) P-wave contribution also features an additional 
(positive) left-right mixing. But most importantly this mixing allows for the smallest 
mass of the f compatible with the lowest mass on the fi set by LEP which is about lOGeV 
lower than for the lightest smuon and selectron. 

First consider the case tan/3 = 10 and r 12 = 0.1. With the constraint Qh 2 < 0.3, one 
finds that at low values of M , the parameter that governs the mass of the right-handed 
sfermions, the lower limit on m^o goes down by a few GeV's relative to the one at large 
values of M , Fig. El However the tail at low M is cut-off to a large extent (by ~ lOGeV) 
when taking into account WMAP. When imposing the constraint flh 2 < 0.3, we find that 
one does not rely exclusively on the Higgsino content of the neutralino to derive a lower 
bound. Large fi values are compatible with neutralinos m^o w 18GeV . However, with the 
much tighter constraint from WMAP, and after including the LEP constraint, the t-channel 
sfermion exchange alone is not sufficient to bring the relic density in the allowed range 
for rrir.o < 3 5 GeV. Once again a non-negligible Higgsino component is required to have 
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pare 3: Region allowed after imposing LEP and relic density constraints for tan j3 = 10, 
= 0.1, [L > and = ITeV, a) in the M — m^,o plane b) in the fi — M Q plane. Here 
scan on M ,M 2 ,/i and A t . Circles (crosses) have flh 2 < 0.3(0.128). Half the points 

the scan are generated in the region corresponding to Mo, /i < 150 GeV. 



sufficient coupling to the Z (Fig. |3J). This entails, as was the case for heavy sleptons, that 
light neutralinos with masses close to the lower bound are necessarily accompanied by 
light charginos. However, for m^o « Mz/2, m-+ can now reach as much as ~ ITeV. 

We show in Fig. H]the lower limit on the neutralino mass for tan (5 = 10 and tan f3 = 50 
as a function of t\i after scanning over M2, n, Mq and A t . The relic density improves 
the constraint on the mass of the light neutralino only after having taken into account 
the LEP2 constraint on selectrons (and on staus for large values of tan/3) as well as 
the constraint from e + e~ — > X1X21 XiXl- After WMAP, the absolute lower limit moves to 
m^o w 29GeV for tan/5 = 10. This limit is obtained in the range r 12 < 0.1 and increases 
with ri2- In the whole region r\i < 0.3, the new WMAP data improves by up to lOGeV the 
lower limit on the neutralino mass obtained with pre-WMAP data. As for heavy sleptons, 
the constraint on the chargino mass from LEP2 sets the lower limit on the neutralino mass 
when 0.3 < ryi < 0.6. We have here also applied the upper bound on the — , Eq. ITHl 
However this additional constraint does not affect the lower bound on the neutralino mass. 
For tan/3 = 50, a lower limit of m^o > 18GeV is obtained for r 12 < .06 and increases 
rapidly with r 12 . 5 This lower limit lies below the one obtained for lower values of tan/3. 
The main reason is the contribution of the channel XiXx ~^ t + t~ with f exchange which 
is enhanced because of the large Tl-Tr mixing. The new WMAP data sets the lower limit on 
neutralino as long as ri 2 < .4. Although we have taken a unified scheme for generating 
the slepton masses along Eq. IHlwith a common Mq at high scale, we would have arrived at 
the same LSP lower limit had we varied the slepton masses independently. This is due, as 
explained earlier, to the fact that from LEP the lowest bound on the slepton applies to the 
f. On the other hand to use the bound from (g — 2) /t to further constrain the parameter 
space tacitly assumes universality, at least as implemented through Eq. El Fig. HJ)) shows 
the impact of applying the (g — 2) M bound alongside the relic constraint. For tan/3 = 10 

5 In [TTJ a lower limit of 18 GeV was obtained for any value of tan/3, using ilh 2 < 0.3, however a 
universal limit on mj > lOOGeV was imposed. 
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Figure 4: Lower bound on the neutralino mass vs r 12 = Mi/M 2 for fi > and a) tan/3 = 
10, b) tan (3 = 50. All direct and precision measurements constraints are combined with 
the LEP limits. We scan over M2,//, Mq and A t . 



there is no impact, however for tan/? = 50 although the lowest bound at very small r 12 is 
hardly affected by (g — 2)^, for intermediate r 12 ~ 0.1 — 0.2 one can improve the limit on 
the LSP bound by about 6GeV. This is related to the fact that both the smallest values 
of Mo as well as of M 2 , corresponding to light smuons and charginos, give too large a 
contribution to (g — 2) M as will be discussed below. 



5.2 fi < 

It has often been claimed that models with negative values of fi are ruled out by the 
(g — 2)^. Within our conservative approach this is not the case not only because the 
present data and the calculation of the hadronic contribution allows <5a /t < but also 
because with nonuniversal gaugino masses the sign of Sa^ and fi are not necessarily cor- 
related. Considering the potential relevance of the (g — 2)^ for /i < models, we will 
first discuss only this constraint together with LEP before folding in the constraints from 
the relic density and b —> 57. To a certain extent we will see that these two sets of 
constraints work in the same direction restricting the region with light sleptons and light 
charginos / neutralinos. 
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5.2.1 (g-2) 



The most important contribution to the (g — 2) M occurs for light smuons and light 
charginos and is enhanced at large values of tan/5. First consider the case r 12 = 0.1 
and light sleptons, M = 150 GeV. When tan (3 = 10, only a few models with very light 
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Figure 5: <5a M in models where \i < 0, r 12 = 0.1 with tan/3 = 10, 50 and a-c) M = 150GeF 
d) M = 1000 GeV. The constraint on 5a^ selects all models above the thick line. Here 
only LEP constraints are imposed. 



large, one derives a lower limit on the neutralino of m^o > 50GeV for tan (3 = 50. It is in 
the small region that one finds the largest deviation in the value of Sa^. For tan (3 = 50 
one quickly exceeds even the conservative 3a bound. However a cancellation between the 
chargino and the neutralino diagrams can change the relative sign of 5a^ and /x, Fig. \^p-c. 
As the neutralino contribution is in general rather small, for this type of cancellation 
to occur, the chargino contribution must be somewhat suppressed (the large M 2 region) 
furthermore a significant amount of Higgsino/gaugino mixing is required, (Fig. Eb(c)). 
Typically this sign flip is compatible with light neutralinos, say < 50GeV, only in the 
context of nonuniversal models where the appropriate hierarchy of parameters can be 
found, that is Mi < \x << M 2 . For heavy sleptons, the predicted values for 5a M are in 
general much smaller and all the parameter space (for r 12 = 0.1) satisfies the experimental 
bound (see Fig. (31). 

Fig.Elshows the exclusion region in the Mq — M 2 plane for three different exclusion val- 
ues for <5a M (depending on the assumptions on the calculation of the hadronic contribution) 
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and for different values of the non universality parameter. Here both \x and A t are kept 
as free parameters and the LEP limits on and on charginos and sleptons are imposed. 
The (g — 2) M rules out regions where M and M 2 (Mi) are small. However, no lower limit 
on the neutralino can be derived as light charginos/neutralinos are allowed when sleptons 
are very heavy. If one would impose the more restrictive limits on £a M , Eq. EE it would be 
possible to accommodate light sleptons and neutralinos in models where r\ 2 < .5, Fig. El 
Indeed, the partial cancellation between the chargino and neutralino diagrams leads to a 
limit on M 2 that is not as strong as in the universal case. With this restrictive bound it 
is possible to set a lower limit on the neutralino from (g — 2) M only for large values of 
tan j3. 
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Figure 6: Contours Sa^ x 10 -10 in the M — M 2 plane with \x < for tan/5 = 10 (left) 
and tan/3 = 50 (right). From top to bottom r 12 = 0.5, 0.1, 0.03. The top panel r 12 = 0.5 is 
based on the usual GUT assumption. Only LEP constraints are folded in. For tan/? = 50, 
the region allowed by b — > 37 lies to the right of the dotted line and as can be seen it does 
not further restrict the bound set by (g — 2) 



iii- 



Note that we have only discussed explicitly the case = 0, while corrections to the 
(g — 2) M are expected for large values of A^ it does not strongly affect the allowed region 
as long as A^ < 1 TeV. Indeed when tan (3 is small there is not much constraint and when 
tan/? is large we have in any case tan/3 >> |/4 M |. 
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Figure 7: Lower bound on the neutralino mass vs r 12 = Mi/M 2 for fi < and a) tan/3 = 
10, b) tan (5 = 50. AZZ direct and precision measurements constraints are combined with 
the relic density limits (dotted lines) while the constraint from Sa^ is removed (full lines). 



5.2.2 Combining all constraints 

The impact of the relic density constraint coupled to direct collider limits from LEP is 
somewhat similar to what was discussed in the case \x > 0. Here we stress the impact of 
adding the indirect constraints from the (g — 2) M . 

• n < 0,tan/3 = 10 

With the new WMAP data, it becomes increasingly difficult to satisfy the relic density 
constraint. Even with light sleptons, sufficient annihilation of neutralinos into fermion- 
antifermion pairs requires the neutralino to be not too far away from the Z pole. As we had 
already found for fi positive the lower limit on the neutralino mass increases significantly 
in the region r%2 < 0.1, from m^o > 18GeV based on the pre-WMAP to m^o > 29GeV with 
WMAP, Fig. [T^i. Note that in this region, the bound on the e + e _ — > XiXi cross section 
from LEP has a significant impact in constraining the small /i region. As a result, the 
lightest neutralino allowed by LEP and WMAP alone occurs for for larger values of the non 
universality parameter, r 12 ~ 0.15 with a bound m^o > 26GeV. It is only by imposing the 
(g—2)^ bound that one increases the lower bound on the LSP in the region 0.1 < rn < 0.4. 
It is also because of the importance of (g — 2)^ that WMAP does not sensibly improve the 
LSP bound in this region compared to the pre-WMAP data. 

• n < 0,tan/3 > 10 

For larger values of tan /3, even our most conservative bound on 5a^ significantly restricts 
the LSP mass, except for the smallest ryi values below about 0.03. Past this value the 
(g — 2) M bound is more restrictive in setting a limit on the LSP mass, see Fig. Eb and 
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very much improves on the WMAP result. As a result the lowest LSP mass m^o ~ 22 GeV is 
found for r 12 ~ 0.03. Had we used the pre-WMAP result, this limit would have been 6GeV 
lower. The effect of (g — 2)^ for r 12 = 0.1 and tan/3 = 50 is more clearly seen in Fig. |HJ 
For this set of parameters m^o > 32GeV are allowed, an increase of nearly 10 GeV had we 
not implemented the (g— 2) M bound. Note that the (g — 2) M bound also cuts on the region 
where both M2(M{) and Mo are smallest. Fig. IHlalso makes it clear that the lightest LSP 
have mass near the Z threshold which as we have argued calls for a not too large \l. In 
turn one sees that this scenario also predicts charginos (with a large Higgsino component) 
that are relatively light, for instance for m^o < 40 GeV, one expects m~+ < 250GeV. 
In the region 0.1 < ri 2 < 0.3, we find that combining all constraints increases the lower 
bound on the neutralino and that the increase is more significant for very large tan/?. 
However, for this range of values for the non universality parameter, WMAP has much less 
impact. In this region the lower limit on the neutralino depends strongly on the lower 
allowed value for 5a^. Had we used the stricter bound on (g — 2) M , Eq. H21 we would not 
have found neutralinos below 60GeV for the large values of tan f3. 
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Figure 8: Allowed region in the M — M 2 plane using direct limits and WMAP constraints 
only (light grey circles) and when taking into account all constraints (dark crosses). Here 
tan j3 = 50, ji < and r i2 = 0.1. 



In summary we have found that the lower limit on the neutralino mass is m^o > 29GeV 
for tan (3 = 10 for either sign of \x and m^o > 18(22)GeV for tan (5 = 50 and /i > 0(/x < 0). 
The latter is allowed only for large mass splitting in the gaugino sector, M\ -C M 2 . The 
muon anomalous magnetic moment does constrain models with \x < as well as the large 
tan/5 region except in the region where Mi -C M 2 . Furthermore, in the region ri 2 > 0.1 
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and when /x < the lower bound on m^o strongly depends on how much negative one 
allows 5a^ to be. 



6 Lowering Ma 



6.1 Limits on the LSP 



The effect of lowering Ma, within the range allowed by LEP, has an impact on both the 
relic density constraint and the indirect limits from the B sector. The relic density can 
be lower than the one obtained in the heavy pseudoscalar case due to the increase of 
the s-channel Higgs exchange contribution (notably the from the pseudoscalar A) to the 
neutralino annihilation cross section. This had to be expected in view of the importance 
of the s-channel pole contribution as we have already seen for the Z and the lightest 
Higgs, h. Recall that in these models Mh and Ma can be as low as 91.6 GeV without 
being in conflict with the LEP2 data. In the large tan (3 region, the contribution of such 
light Higgses can be sufficient to bring the relic density below the upper bound even if 
one is quite far from the s-channel pole due to the enhanced coupling of the Higgs to 
b quarks and r's. This opens up a new region of parameter space for light neutralinos 
(m^o < 18GeV) as was pointed out in [73 EH]- 




Figure 9: Vth 2 as a function of M A for fj, = 120GeV, M = 1000 CeV, M x = 0.01M 2 
and a) tan/5 = 10 b) tan/5 = 50. Full lines correspond to (from top to bottom) M 2 = 
800,1200,1800 GeV (left) and M 2 = 300,500,800,1200,1800 GeV (right). Additional 
dashed lines show the effect of varying /j when M 2 = 800 GeV (from bottom to top) // = 
160, 200, 240, 280Ge V. The dash line indicates the current upper bound from WMAP. Only 
direct LEP limits are imposed. 



The dependence of the relic density, for small /j, on Ma is displayed in Fig. |H1 for 
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r 12 = 0.01 and for tan/3 = 10 and tan/3 = 50. The upper bound on the relic density is 
satisfied most easily for a light pseudoscalar, Ma ~ 100 — 120GeV and for large values 
of tan/3. For instance tan/3 = 10 is totally excluded by WMAP for the set of parameters 
chosen, in particular for heavy sleptons, as seen in Fig. Ek- Due to the enhanced coupling 
of the pseudoscalar to the heavy fermions, oc tan/3, Fig. E3 also shows that changing tan/3 
from 10 to 50 makes the relic density compatible for low values of Ma- The relic density 
decreases as one increases M\ (or Mi since rxi is fixed), thus increasing the LSP mass. 
Once again, we confirm that in order to maximize the Higgsino-bino mixing, small values 
of fi are preferred. Generally, the light neutralinos (say below 16GeV) are found when 
Mi « n « M 2 and we therefore also expects relatively light charginos as was the case 
with the heavy pseudoscalar scenario. 

In models with low Ma and large tan (3 one expects large effects from the b observables. 
Most affected are b — > s'j and B s — > fi + fi~ . On the other hand we have checked that the 
branching ratio to Z — > bb remains insensitive even in this scenario. The branching ratio 
for b — > 57 receives tan f3 enhanced contributions from both the charged Higgs diagrams 
as well as the chargino/stop diagrams. Typically the contributions form each types of 
diagrams individually could be far above the allowed limit. However a partial cancellation 
between these two contributions occurs when A t is large and negative. For example for 
tan/3 = 50, one needs values of A t < —750 GeV. For tan/3 = 10, on the other hand, even 
the value A t = —2400 GeV, is not sufficient to induce a cancellation between the chargino 
sector and the large contribution from the charged Higgs to b — > S7. This process then 
forbids the very low values of Ma for intermediate values of tan /?, Fig. One sees for 




m x o(GeV) 

Figure 10: Allowed values for m^o as a function of the pseudoscalar mass for tan (3 = 10. 
The impact of the WMAP data on the relic density is displayed: pink (light grey) circles 
corresponds to Qh 2 < 0.3 and green(medium grey) crosses to Qh 2 < 0.128. The black 
crosses correspoind to the region allowed after including the b — > 37 constraint. 

example that the lower bound for tan/3 = 10, moves from m^o pa 12GeV to 29GeV after 
imposing the b — > 57 constraint. Models with very light neutralinos that are compatible 
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with the .B-sector constraints as well as with the relic density are then expected to be 
models with light pseudoscalars and large values of tan (3 as one can see in Fig. [TTJ 
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Figure 11: Region of the taxi (3 — Ma plane where one can find light neutralinos, Light grey 
(green) crosses have 10GeV< m^o < 16 GeV while black crosses correspond to m^o < 
10 GeV. Here all constraints are implemented. 



In these models, the branching ratio for B s — > /x + can be strongly enhanced at low 
values of Ma- In particular diagrams with chargino loops could give a large contribution 
when the chargino has a large Higgsino component. This corresponds to not too large 
values of M 2 . The predicted values for B s — > fi + fi~ are displayed in Fig. ^] for a variety 
of models that pass both LEP and WMAP constraints. One sees that the present limit from 
the Tevatron on the B s — > / u + [i~ eliminates some models but leaves open the possibility 
of light pseudoscalars. 

A scan over the full parameter space, Eq. with 92GeV < Ma < 500GeV and 
imposing all the above mentioned constraints leads the lowest LSP mass bound of m^o = 
6 GeV. 6 This LSP is found in models with large values of tan f3 tan (3 > 30 and is 
associated with Mh, Ma, \fJ>\ < 120GeV. The sleptons need not be light. The allowed 
region in the Ma — rn^o plane (Fig. IT3*|) for tan (3 = 50 clearly shows that when Ma < 
250 GeV one can lower the bound on the neutralino compared to the case when the 
pseudoscalar is 1 TeV. Contrary to what we have seen for other cases, the impact of the 
WMAP results is marginal as concerns the lower bound on the neutralino mass ~ 6 GeV, 
see Fig. Although we found that a large number of models were ruled out by the 
requirement of b — > S7 fEq. ITTj). this constraint does not affect the absolute lower limit on 
the neutralino mass, m^o = 6 GeV, in the large tan (3 scenario. Similarly the B s — > fi + fi~ 
constraint does not impact on this lower bound. Moreover, for intermediate values of 
tan f3 ~ 10 — 20, we find that lowering Ma does not affect the lower limit on the neutralino 
mass which we derived for Ma = 1 TeV. This is due to an incompatibility with the b — > S7 
constraint as shown explicitly for tan/5 = 10 in Fig. El In fact, restricting the analysis 
to the region where m^o < 16GeV, we found that only models with tan/3 > 25, \fi\ < 

6 This agrees with the results of |88j . 
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Figure 12: Branching fraction for B s — > fi + /j,~ vs the pseudoscalar mass for tan/3 = 50, 
\i > 0. We scan over all parameters and impose the LEP and WMAP constraints (green/light 
grey) as well as all constraints (dark). The preliminary 95% upper limit obtained by 
CDF with 113 pb^ 1 (9.5 x 10~ 7 ) and the expected reach of Runlla (2 x 10~ 7 ) are also 
display edJTSlj. 
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Figure 13: Allowed values for m^o as a function of the pseudoscalar mass for tan (3 = 50. 
The impact of the data on the relic density is displayed. Green stars (medium grey) 
correspond to Qh 2 < 0.3 while black crosses correspond to Qh 2 < 0.128. b) Chargino 
masses as a function of LSP mass. All constraints are implemented. 
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400GeV, Ma < 200GeV and a large negative stop mixing parameter A t < —750 GeV 
were consistent with all constraints. The first three conditions are necessary for sufficient 
annihilation of the LSP into rf(bb) pairs, while the last condition ensures cancellation 
between the chargino and charged Higgs contributions to b — > sj. Obviously, as the 
coupling of the light Higgs scalar is enhanced with tan/3, the range of Higgs masses 
allowed is wider for larger values of tan/3, (see Fig. [TTJ). With not so heavy pseudoscalar 
Higgs and charginos, colliders, and in particular the Tevatron, might then have good 
prospects to constrain the models or discover new particles in such a framework as will 
be discussed next. 



6.2 Prospects at the Tevatron 

Based on the existing Tevatron analyses |74j. we infer that the Tevatron should be sen- 
sitive to some neutral Higgs bosons in a large fraction of the models having very light 
neutralinos with m^o < 16 GeV. As we have just discussed, the various constraints already 
impose that Ma < 220 GeV and that tan /3 be large to arrive at such low values for the 
neutralino mass (see Fig. ITT]) . Then the channel pp — > bb<f), with — > bb for <p = h, H, A 
is rather sensitive to the large tan /3 region due to the enhanced coupling to the b quarks. 
For example, for tan/3 = 50, in the Mh max scenario, values of Ma ~ 105GeV are excluded 
already with 0.1/&~ 1 of integrated luminosity and with 2/6 -1 the region Ma < 200GeV 
can be excluded at 3a [74]. For tan/3 = 40(30), M A < 160(130) GeV can be excluded 
with 2fb~ l . For the very large values of tan/3 under consideration we have checked that 
the invisible mode, h — > XiXi carries only a small fraction of the total width so that the 
bb branching ratio would not be significantly suppressed. The models with light Higgses, 
Mh,H,A ~ 100GeV can also be probed at the Tevatron in the charged Higgs channel 
via disappearance searches for light charged Higgses produced in the decay of top quarks 
t — > H + b. With 2fb~ 1 of integrated luminosity in Runlla, the region m#+ < 130(150)GeV 
for tan/3 = 30(50) can be excluded at 95%. However, as the charged Higgs are heavier 
than the pseudoscalar, this channel does not improve on the potential of the Tevatron . 

Moreover the required values of \l for models leading to m^o < 16GeV are such that 
the lightest chargino masses can also fall within the range accessible by the Tevatron in 
Runlla. In particular, when Ma gets close to 200GeV and higher luminosities are needed 
for the Higgs channels, the predicted values for the charginos cluster near 100 — 150GeV. 
Although a detailed analysis of chargino searches in the trilepton mode in non-universal 
models has not been completed yet, [34], it is expected that the Tevatron would be 
sensitive to charginos within this range. 

Finally, the Tevatron can also probe the large tan /3-light LSP scenario, via the B s — > 
fi + pT . In Fig. H2J we show the range of predicted values for B s — > as well as the 

expected sensitivity of Runlla (2 x 10~ 7 ). A large number of models with light neutralinos, 
predict a branching ratio above this sensitivity. The largest branching ratios are found 
for M A < 250GeV. 

To estimate the potential of Runlla to probe models with very light neutralino LSP, 
we assume the expected limits, rriA < 130(200) GeV, m H + < 130(150) GeV and B s — > 
fi + {i~ < 2 x 10~ 7 , for values of tan/3 = 30(50). For the chargino channel we take as 
a guideline the value m-+ < 150 GeV. Basically, we find that the additional region of 
parameter space where light neutralinos were allowed as one lowered Ma will be probed 
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entirely, at tan/3 = 50, by the Tevatron , both in the pseudoscalar searches as well as in 
the B s — ► The chargino searches will also somewhat help constraint the allowed 

models. For tan/3 = 30, due to the lower reach in the Higgs channel, we found some 
models with m^o > 15GeV which could escape detection in Runlla. 

In summary, decreasing the value of M4 opens up new possibilities for light neutralinos 
(6 — 16GeV). With Higgs and chargino searches, the Tevatron should be able in the near 
future to cover the small remaining part of parameter space where very light neutralinos 
(~ 6GeV) can exist. 

6.3 Direct detection with light scalars 

The models we have just discussed admitting pseudoscalars (and hence the other neutral 
Higgses) with masses close to the lowest limit set by LEP offer interesting prospects for 
the direct detection experiments as well. The spin independent scalar cross section for 
neutralino scattering on a nucleus is sensitive to squark and Higgs exchanges. However, 
in the models with heavy squarks considered here, the dominant contribution arises from 
t-channel Higgs exchange. The dependence of the cross section on the mass of the Higgs 
scalars goes simply as a^- oc j^r. For large values of tan/3, the enhanced coupling of the 

Higgs to d-type quarks tends to further increase the cross section. Moreover because of 
the low /i values that these models admit, the LSP coupling to the Higgses is maximised. 
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Figure 14: The spin independent neutralino-proton cross section as a function of a) Ma 
and b) m^o for tan/3 = 10 (dark grey) and tan/3 = 50 (green/light grey). Scan over 
the parameters specified in Eq. All constraints were imposed. The reach of the di- 
rect detection experiments ZeplinI (full) , Edelweiss (dot) CDMS (dot- dash) as well as for 
ZeplinI I (dash) is displayed. 

We have calculated the scalar cross section for neutralino scattering off nucleons. We 
have included the contribution of all six flavors of quarks. Our results agree with the 
ones given in Ref . [T^] . We have used the values for the coefficients of the quark mass 
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operator in the proton, frf* as given in Ref. [ZH]- We show the predictions for the spin- 
independent cross section for neutralino-proton scattering as a function of Ma, Fig. Hlk . 
after scanning over the parameter space as specified in Eq. EH One witnesses the large 
enhancement at low values of Ma- The various constraints from LEP, cosmology as well 
as from precision measurements have been folded in the global scan over the parameter 
space, hence the much smaller predictions for the cross section at intermediate values of 
tan (3. There, as was discussed above, the b — > 57 constraint is not compatible with light 
pseudoscalars. The cross section for tan f3 = 10 therefore never exceeds a few 10 _8 pb. 

For tan/? = 50, one reaches cross sections that are in the region of detectability of 
CDMS/Edelweiss/Zeplin although the largest cross sections are often found in regions 
where neutralinos are very light, precisely where detectors have shown less sensitivity, 
Fig. 114b . Nevertheless a few models, albeit not the ones with the LSP near the lowest 
limit ps 6GeV, are already ruled out by CDMS and Zeplin[29j. Furthermore, in the majority 
of models with M4 < 200GeV, the spin-independent cross section exceeds 10~ 8 , detectable 
by the next generation of direct dark matter detectors such as EdelweissII or Zeplinll. 7 
We emphasize that to cover all models it is crucial to have a good sensitivity in the region 
m^o « 6 — 10 GeV and that future detectors should extend their capabilities and searches 
in this mass range. 




H (Gev) 

Figure 15: Regions of the Ma~^ plane where the spin independent neutralino-proton cross 
section exceeds 10~ 7 pb (black), 10~ 8 pb (green/light grey), 10~ 9 pb (blue/medium grey). 
Here, tan (3 = 50 and we scanned over the parameters specified in Eq. 03 All constraints 
were imposed. 



7 Note that there are uncertainties in the exclusion curves written in terms of the LSP scattering cross 
section on proton, for a discussion of the effect of halo modeling and LSP velocity distribution see|77|. 
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Although the spin independent cross sections are enhanced at large tan/?, cross sec- 
tions above 10 _7 pb are found for all values of tan/3 where light pseudoscalars are allowed 
(tan/3 > 25). Isocurves of the scalar cross section are displayed in the Ma — p plane 
in Fig. emphasizing the importance of a sizable amount of gaugino-Higgsino mixing 
to obtain large spin-independent cross sections. Here only the upper bound on the relic 
density has been considered. Many of the models, near the Z or Higgs resonance actually 
have a value for the relic density that is much too low. Usually one rescales the cross 
section by a factor of the relic density to some minimal value(say 0.05) in order to take 
into account the fact that the neutralino would not constitute the main cold dark matter 
component. Applying a rescaling factor of p = Qh 2 /0.05, would strongly affect the effec- 
tive cross section(<j e // = pc^p') near * ne % or Higgs pole but not for the very light LSP 
where the relic density is near the allowed upper bound. 



7 Prospects at e + e colliders 

The aim of this section is to address the issue whether SUSY models with non-universal 
gaugino masses that have a light LSP, say m^o < 70GeV, can always be probed at an e + e~ 
linear collider of centre of mass 500GeV. In the usual mSUGRA type model, imposing 
such a light neutralino entails m~+ < 150 — 200 GeV and therefore independently of any 
other sparticles, chargino discovery is guaranteed [THj- The non-universal models do not 
always ensure m~+ < 250 GeV, moreover some models do not have light sleptons as we 
have seen and yet are consistent with all data (including cosmology). 

The main channels for producing supersymmetric particles at a 500GeV linear collider 
are the slepton pair production and chargino pair production with e + e~ —>■ XiXj being 
important processes as well. Of course the lightest Higgs boson can always be produced. 
In non universal models, the difficult cases for the LC will be the ones where sleptons 
are heavy with the mass of the LSP near Mz/2 S . In this case large values of p are still 
allowed. This entails charginos (and X2) t°° heavy to be directly pair produced. The only 
processes accessible would then be XiX,] 

To ascertain which (s)particles can be reached at a 500 GeV, we perform scans over 
the parameters M 2 , p, M for fixed values of tan f3 and r 12 , the non universality parameter. 
We compute the unpolarized cross sections for e + e~ — > tRtR, I^bP-r-, T\f\, XiXi > X1X2; X1X3 
as well as for XiXil- We have not simulated ex,e# and e^e^ production. In our scenario 
these would be accessible only when e^e# is accessible. The cross sections are computed 
with GalcHEP|41[ I79j. a program for the automatic calculation of Feynman diagrams in 
the standard model or the MSSM. Because the LSP is very light, clean signatures for all 
sparticles that decay into Xi are guaranteed. In fact, there is in general sufficient phase 
space for the decay into LSP to constitute the main decay mode for the selectron, — > 
e~Xi, the chargino xt ~~ * ^ + Xi as wen as heavy neutralinos, Xi ~~ *" Xi^^- F° r charginos 
and neutralinos, decays into sleptons are also sometimes kinematically accesssible. We 
will not perform a detailed analysis of signal and background for all these processes, in 
the e + e~ clean environment it should not be a problem to see a signal for a > lib. 

8 These models should not be thought of as fined tuned, from the relic density point of view, more 
than the mSUGRA model is. The latter for example does require, for example, near degeneracy between 
the fi and Xi masses to pass the WMAP constraint. 
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We also consider e + e~ — > XiXili the process that contribute to a 1+inv . Although this 
process has to fight a large background from the standard e + e~ — * uu^y, with 95% e~ 
polarization and after a cut on the small angle photons, it was shown that it was possible 
to reduce the background to < 150fb.[BD] From this we estimate that with a luminosity 
of 500/6 -1 , the polarized cross section for the signal should exceed a = 1.6 fb in order 
to see a 3a deviation. For simplicity we compute only the unpolarised cross section and 
impose conservatively the same value as the pair production processes, cr unpo i > lib. 

Here we concentrate on the \x > region. We find that in large regions of the parameter 
space the cross section for at least one process exceeds lib. We do, however, uncover 
regions compatible with existing limits where no sparticle is produced. To look more 
precisely at the potential of a 500GeV LC, we discuss again two typical values of tan j3. 

• tan/5 = 10 

First consider the case tan/3 = 10. For unified models (r 12 ~ .5) and with m^o < 
70 GeV, the whole parameter space is covered basically by the processes e + e~ — > XiXi ■ 
While chargino pair production is the preferred channel especially for the large Mq region, 
in this scenario, X2X1 a l so exceeds lib in the full parameter space, Fig. [TBI The process 
X3X1 has a significant cross section only for // < 400GeV. When M < 200GeV, the 
slepton pair production is also always accessible. As one decreases the value of the non 
universality parameter, r 12 = 0.1, one starts uncovering permitted regions where none of 
the charginos/neutralinos production processes have a large enough cross section while 
the sleptons are too heavy for pair production. This occurs in the large Mq region and for 
very large \x where the coupling to the Z is not sufficient. As the chargino mass can be 
much higher than the LSP, the chargino process is accessible only in the small \i region 
while the X1X2 an d X1X3 both extend the reach in parameter space. The process X1X2 * s 
also observable in a narrow region at small M even for large values of these were the 
points where m^o ~ Mz/2. 

When ri2 decreases to 7*12 = 0.03, one observes the same features. For M < 200GeV, 
observable cross sections for the slepton production processes crCr, [Ir^r and fx fx are 
expected as for the universal case. For larger values of M , it is o^o^o that eventually 
exceeds cr^o^o in the region M\ « fi < M 2 . Then the LSP is bino-like and both X2 an d X3 
are Higgsino-like and nearly degenerate. However the ZxiXs features the largest coupling. 
With the new data from WMAP, the large Mo/// region is more severely constrained, thus 
reducing the allowed parameter space where no supersymmetric particles can be produced. 

Note that one can also get indirect evidence of light neutralinos by measuring the 
invisible branching fraction of the Higgs. A large branching fraction into invisible is ex- 
pected in models with non-universal gaugino masses in the region of intermediate tan/3 
and when /x is small p33]- For example, for tan /3 = 10, n 2 = 0.1, /i < 300 GeV the branch- 
ing fraction h — > X1X1 can reach 60%. This region of parameter space corresponds to the 
one where one expects many other signals of supersymmetry at the LC500, XiXi>XiXj 
and sometimes crCr, Fig. Correlating the information from the Higgs and the super- 
symmetric sector will then help to establish the nature of the invisible decay mode of the 
Higgs. 

• tan/? = 50 

In the large tan/3 scenario, the slepton pair production process is not so useful, at 
least in the unified model, as most of the low M region is ruled out. The selectron pair 
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Figure 16: Region of Mq-/i plane where the cross section for sparticles production in 
e + e~ exceeds Ifb, for tan/5 = 10 (left) and tan/? = 50 (right) and from top to bottom 
r 12 = 0.5,0.1,0.03. r 12 = 0.5 is the usual GUT value. The processes are crCr (yellow 
hatch) XiXi (red/medium grey) X1X2 (blue/dark grey) X1X3 (green/light grey). Note that 
there is a large overlap between the region covered by crCr and the ones (not featured) 
covered by fiRpiR or f\f\. The region to the right of the full (dashed) line do not satisfy 
the constraint Qh 2 < 0.128(0.3) coupled with LEP and precision measurements. In the 
region (white) to the left of the WMAP line no supersymmetric particle can be produced. 
Such regions are present only when r\2 = 0.03 and ryi = 0.1. In the top figures, the whole 
parameter space is also covered by XiX.2- ^ n ^ e bottom right figure, the dotted line delimits 
the region where o~(e + e~ — > XiXil) > 1/^- 



29 



production process is only accessible in a corner of parameter space when /j,, M m lOOGeV. 
The stau pair production is also accessible when \i is large enough to induce significant 
mixing in the f sector. However, with the WMAP constraint there remains only a very 
narrow region at large \x. Both chargino pairs and X1X2 cover all allowed parameter space 
while X1X3 is al so important when \x < 400 GeV. As rvi decreases, the low Mo region is 
allowed and observable cross sections for tRtR as well as f\f\ are predicted. Again there 
is a portion of the large Mo-large \i region of the parameter space where supersymmetric 
particles are too heavy to be produced. Note that imposing the new WMAP constraint 
has considerably shrunk this region. Finally, the cross section for the process X1X1I can 
exceed lib, although this occurs only in the low M region when slepton pair production 
is also accessible. 

Increasing the energy of the collider would obviously increase the parameter space 
where sparticles could be pair produced. Basically, the slepton pair production is large 
for M < 400GeV and chargino pair production for // < 450GeV. Thus, one nearly 
recovers full coverage with only these two processes, in the case of nonuniversal models 
with light neutralinos at a 800GeV collider. At first sight, the more problematic models 
are those with ryi ~ 0.1 where much larger values of \i are allowed at intermediate tan/3. 
However, the process e + e~ — > X1X2 nas a reasonable cross section for most of the large 
Mo — /i region. For example when tan/3 = 10, X1X2 i s measurable everywhere except in 
a small region where Mo > 900GeV. Then both the t-channel sfermion exchange and the 
s-channel Z exchange contributions are small and no signal for supersymmetric particles 
can be expected. With WMAP, this difficult region is already ruled out when r 12 < 0.1 
(see Fig. El for the allowed region). One also avoids the difficult region when tan/3 = 50 
(Fig. ITfij) . There, for any value of the non universality parameter the allowed region 
corresponds roughly to the one where chargino pair production is accessible. 

Finally we mention that such light supersymmetric particles would probably first be 
discovered at the hadron colliders. For the Tevatron , although the chargino mass con- 
sidered are often within the range where discovery via the trilepton signal is expected 
(m-+ up to 250 GeV) one needs to take a closer look at cross sections for signal and 
backgrounds in nonunified models to ascertain the viability of the signal [33]. In the 
models with heavy squarks and Higgs scalars, no other opportunity for SUSY discovery 
exists. The LHC on the other hand would have plenty of opportunities to discover SUSY 
in either the squark, Higgs or gaugino sector. We should also note that the LC through 
precision measurements on the h properties could also probe into SUSY especially by 
using information from the LHC. 

8 Conclusion 

We have reexamined the lowest bound on the neutralino LSP in the minimal supersym- 
metric model. We have worked within the context of minimal flavour violation and R 
parity conserving supersymmetric models. We have reduced the parameter space to only a 
few important parameters, those of the gaugino and slepton sector, the pseudoscalar mass 
and the trilinear coupling of the squarks. In particular we have relaxed the universality 
relation between the gaugino masses thus removing the most important constraint on the 
LSP mass arising from LEP. We find that the upper limit on the relic density contributed 
by neutralinos and as inferred from the new data by WMAP basically sets the lower bound 
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on neutralinos in models where the gaugino masses are not unified at the high scale but 
satisfy Mi << M 2 . In the limit of a heavy pseudoscalar mass, we found a lower bound on 
the neutralino mass of 18GeV in models with tan /3 = 50. For intermediate values of tan /3, 
the annihilation of neutralinos into light fermions is not as efficient and one can only set 
a bound of 29GeV for tan (3 = 10. It is however in models with light pseudoscalar masses 
(Ma < 200GeV) and large values of tan/3 that one finds the lightest LSP m^o « 6GeV. 
This is due to a new contribution to the annihilation cross-section of neutralinos, scalar 
exchange mediating XiXi t + t~ , thus reducing the relic density. 

The models with light pseudoscalars and neutralinos below « 16 GeV can be probed 
both in the next generation of direct detection experiments as well as at the Tevatron or 
future colliders. The cross section for the former are much enhanced for a light Higgs with 
significant coupling to the neutralino. It is however crucial that the detector be sensitive 
to light neutralino masses. In addition, the pseudoscalar and/or scalar Higgses should be 
within reach of the Tevatron Runlla. Charginos or B s — * might also be detected 

there. The LHC also has a potential for discovering charged Higgses in the large tan/3 
region using gb -> tE ± with H + -> tb, t + u or qq — > H ± — > rv |81j . Studies of neutral 
Higgses searches in the t + t~ channel, also show a good potential for discovery even with 
a low luminosity |81j . Of course other sparticles, such as charginos, heavier neutralinos 
or sfermions might be discovered as well. 

Finally we have shown that a linear collider with centre of mass of 500GeV has good 
prospects for producing supersymmetric particles in models where neutralinos are below 
the weak scale, although the parameter space cannot be completely covered especially in 
models where M\ << M%- The gaugino universality relation can be directly tested at such 
a collider but this necessitates performing a combined fit to the neutralino mass (measured 
in slepton pair production process), the chargino mass as well as to the polarized cross 
sections for e + e r R — > e^ejj, XiXi jHUlEI] All these precision measurements can be realized 
only if M (m gl ),/x < 250GeV. Of course models where Ma is near (9(100) GeV should 
lead to several signals at the linear collider. Most noticeably in the Higgs sector where 
one expects a large cross section (10 — 100/6) in both e + e~ — *> ZH, hA channels as well 
as in the charged Higgs pair production channel. In most cases, and most importantly 
when the Higgs sector is harder to probe, the charginos would be accessible as well at a 
collider with <fe = 500GeV. 
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